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A  micromodel  based  upon  percolation  theory  is  developed  to  predict  effective  properties  in  composite 
electrodes  for  solid  oxide  fuel-cell  (SOFC)  applications.  The  theory  considers  binary  and  multi-component 
mixtures  of  particles  that  are  either  ion  or  electron  conductors.  The  model  predicts  effective  ionic  and 
electronic  conductivities,  three-phase  boundary  lengths,  and  hydraulic  pore  radii.  The  effective  proper¬ 
ties  depend  upon  primary  physical  characteristics,  including  average  particle-radii,  volumetric  packing 
densities,  particle  contact  angles,  and  porosity.  All  results  are  presented  in  nondimensional  form,  which 
provides  considerable  generality  in  their  practical  application. 

©  2009  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

composite  electrodes  play  critical  roles  in  most  solid  oxide 
fuel-cell  (SOFC)  architectures.  A  typical  composite  anode  is  a  mix¬ 
ture  of  Ni  and  yttria-stabilized  zirconia  (YSZ).  A  typical  composite 
cathode  is  a  mixture  of  lanthanum-doped  strontium  manganate 
(LSM)  and  YSZ.  It  is  well  known  that  the  composite-electrode 
architecture,  including  functional  grading,  is  important  to  overall 
fuel-cell  performance  [1-3].  Recently,  several  groups  have  reported 
significantly  improved  performance  by  introducing  nanoscale 
particles  into  composite  electrodes  [4-8].  There  is  certainly  a 
need  to  understand  quantitatively  the  influence  of  composite 
microstructure  on  fuel-cell  operation.  Applying  that  understanding 
is  important  to  the  design,  optimization,  and  fabrication  of  certain 
microstructure  architectures  as  a  means  to  achieve  high  fuel-cell 
performance. 

Composite  electrodes  are  frequently  thought  of  and  mod¬ 
eled  as  random  mixtures  of  spherical  particles.  The  mixtures 
contain  electrolyte  particles  (e.g.,  ion  conducting  YSZ)  and  elec¬ 
trode  particles  (e.g.,  electron  conducting  Ni).  Composite-electrode 
fabrication  usually  begins  by  mixing  together  electrode  and 
electrolyte  particles  followed  by  sintering  at  high  temperature. 
High-resolution  microscopic  images  of  actual  composite  electrodes 
reveal  that  the  “particles”  are  not  spherical,  even  when  the  particles 
may  have  been  initially  spherical  before  sintering  [9].  Neverthe¬ 
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less,  models  based  upon  random  arrays  of  spherical  particles 
provide  important  qualitative  insight  and  quantitative  design  guid¬ 
ance. 

The  present  paper  develops  and  extends  coordination-number 
and  percolation  theory  to  predict  effective  properties  of  compos¬ 
ite  electrodes.  Based  upon  theory  at  the  particle  scale,  important 
characteristics  of  the  composite  structure  can  be  derived.  Such 
effective  properties  (e.g.,  electron  and  ion  conductivity,  three- 
phase  boundary  length,  and  pore  radius)  are  needed  in  larger  scale 
models  that  consider  the  entire  fuel-cell  system  [10,11  ].  The  effec¬ 
tive  properties  depend  upon  microstructural  features,  including 
particle-radii  distributions,  gas-phase  porosity,  particle  overlap, 
and  relative  electrode  and  electrolyte  particle  loading.  The  com¬ 
posite  performance  also  depends  upon  intrinsic  properties  (e.g., 
conductivity)  of  the  particle  materials.  The  results  herein  are 
presented  in  nondimensional  form,  affording  a  degree  of  general¬ 
ity. 

The  percolation  theory  in  this  paper  builds  on  the  founda¬ 
tion  of  significant  earlier  literature.  As  part  of  a  more  general 
effort  to  develop  micromodels  of  SOFC  systems,  Costamanga  et 
al.  [12,13]  developed  percolation  models  for  composite  electrodes. 
Chan  et  al.  [14,15]  also  developed  SOFC  micromodels  that  incorpo¬ 
rate  percolation-theory-based  effective  properties.  Janardhanan  et 
al.  [16]  used  percolation  theory  to  predict  three-phase  boundary 
lengths  as  functions  of  particle  dimensions  and  packing  den¬ 
sity. 

Percolation  theory  depends  upon  the  concept  of  coordination 
numbers,  representing  the  number  of  contacts  that  a  certain  parti¬ 
cle  makes  with  neighboring  particles.  The  micromodel  developed 
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Nomenclature 

ak,l 

cross-sectional  area  per  contact  between  a  k-  and  an 
^-particle  (m2) 

aV 

ed.el 

overall  binary  phase  contact  area  per  unit  volume 

(m-1) 

interfacial  area  per  unit  volume  between  the  solid 
particles  and  pores  (nrr1 ) 

ed  k 

subscript  for  electrode  particles  in  size  class  k 

elk 

subscript  for  electrolyte  particles  in  size  class  k 

ed 

subscript  for  m  kinds  of  electrode  particles 

el 

subscript  for  n  kinds  of  electrolyte  particles 

contact  perimeter  between  a  /c-particle  and  an  t- 
particle  (m) 

L 

geometric  thickness  of  composite  electrode  (m) 

Lf 

total  thickness  of  grain  boundary  along  the  current 
path  (m) 

M 

number  of  particle  types 

m 

kinds  of  size  class  for  electrode  particles 

< 

number  of  /c-particles  per  unit  volume 

ni 

number  of  /c-particles  per  unit  electrolyte  surface 
area 

n 

kinds  of  size  class  for  electrolyte  particles 

Nk,i 

average  contact  number  of  ^-particles  with  a  k- 
particle  in  an  idealized  medium 

h 

percolation  probability  of  /c-particle 

rk 

radii  of  /c-particles  (m) 

rc 

neck  radius  between  two  connect  particles  (m) 

rg 

hydraulic  radius  for  gas  flow  (m) 

h 

nondimensional  hydraulic  radius 

Re|r 

ion- transport  resistance  (£2) 

h 

surface-area  fraction  of  all  ^-particles 

s 

geometric  cross-sectional  area  of  composite  elec¬ 
trode  (m2) 

sf 

effective  contact  surface  area  between  el-particles 
per  particle  layer  (m2) 

X 

fraction  of  solid  surface  areas  belonging  to  the  solid- 
pore  binary  phase  boundary  area 

zk 

average  number  of  contacts  between  a  /c-particle 
and  neighboring  particles  of  all  types 

Zk,i 

average  number  of  contacts  between  a  /c-particle 
and  ^-particles 

Z 

overall  average  coordination  number  of  all  solid  par¬ 
ticles 

Ze\,e\ 

average  coordination  number  for  all  n  electrolyte 
particles 

Greek  letters 

8 

thickness  of  boundary  between  particles  (m) 

€ 

relative  error  between  ^kZk  i  and  ^Zik 

\V 

ATPB,eff 

TPB  length  per  unit  volume  (m2) 

effective  TPB  length  per  unit  volume  (m2) 

ATPB,eff 

nondimensional  effective  TPB  length  per  unit  vol¬ 
ume 

}S 

ATPB,eff 

effective  TPB  length  per  unit  electrolyte  surface  area, 

rr»-l 

ss 

ATPB,eff 

m 

nondimensional  effective  TPB  length  per  unit  elec¬ 
trolyte  surface  area 

fl 

Bruggeman  factor 

1* k 

volume  fraction  of  /c-particles  relative  to  the  total 
solids 

n 

percolation  threshold  of  /c-particles 

redk 

relative  volume  fraction  of  /c-c lass  particles  within 
all  m  electrode  particles 

4>g 

porosity 

CTtra’° 

intrinsic  material  conductivity  (S/m) 

^tra.eff 

effective  intra-particle  conductivity  (S/m) 

^jtra.eff 

nondimensional  effective  intra-particle  conductiv¬ 
ity 

(jter.o 

intrinsic  conductivity  of  the  particle  interface  (S/m) 

^ter.eff 

effective  inter-particle  conductivity  (S/m) 

^jter.eff 

nondimensional  effective  inter-particle  conductiv¬ 
ity 

0 

contact  angle  between  two  particles 

£el 

effective  relative  density  of  electrolyte  particles 

(k 

number  fraction  of  /c-particles 

by  Bouvard  and  Lange  has  been  widely  used  to  estimate  coordina¬ 
tion  numbers  for  binary  systems  [17].  Suzuki  and  Oshima  developed 
a  theory  that  can  be  used  to  estimate  coordination  numbers  in 
multi-component  mixtures  [18-20]. 

The  present  paper  develops  some  significant  extensions  of  ear¬ 
lier  research.  The  approach  considers  multi-component  particle 
mixtures  and  satisfies  the  contact-number  conservation  require¬ 
ment  [12].  Based  upon  the  underpinning  functional  forms,  all  the 
effective  properties  are  represented  in  nondimensional  form.  Con¬ 
sequently,  the  results  are  general  in  the  sense  that  they  are  not 
specific  to  a  particular  electrode  structure.  Rather,  once  physical 
parameters  are  specified,  relevant  effective  properties  can  be  easily 
extracted  from  the  nondimensional  properties. 

It  should  be  mentioned  that  direct  simulation  of  large  particle 
arrays  is  an  alternative  to  percolation  theory  for  estimating  effective 
properties  [21-25].  Additionally,  such  multi-particle  simulations 
can  provide  quantitative  information  concerning  the  coordination 
numbers  that  appear  in  the  percolation  theory. 

2.  Percolation  theory 

Percolation  probabilities  represent  the  likelihood  that  particles 
are  clustered  in  ways  that  form  connected  conduction  pathways, 
which  is  central  to  the  effective  functioning  of  a  composite  elec¬ 
trode  [3,13,21].  The  concept  of  coordination  numbers  is  central 
to  the  practical  implementation  of  percolation  theory.  Broadly 
speaking,  the  coordination  number  represents  the  number  of 
contacts  a  particular  particle  makes  with  its  neighboring  parti¬ 
cles. 

Consider  a  multi-component  mixture  of  randomly  packed 
spherical  particles  with  M  particle  types  (M  is  typically  a  small  num¬ 
ber).  For  example,  in  a  binary  mixture  of  Ni  and  YSZ  particles  each 
with  a  single  particle  radius,  M  =  2.  The  average  number  of  contacts 
Zk  between  fc-type  particles  and  neighboring  particles  of  all  types 
is 

M 

zk  =  )Tzm,  (1) 

1= 1 

where  Zki  is  the  number  of  contacts  between  a  /c-particle  and  t- 
particles.  Suzuki  and  Oshima  proposed  that  Zki  is  proportional  to 
Si  (the  surface-area  fraction  of  all  ^-particles)  and  Nk  i  (the  average 
contact  number  of  all  ^-particles  with  a  /c-particle  in  a  particular 
idealized  medium)  [18,20]: 


zk,l  =  SiNk,i. 


(2) 
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The  surface-area  fractions  can  be  estimated  in  terms  of  the  particle 
volume  fractions  V'k  as 

n\Anr\  fcrf 

S ‘  =  ~u - =  ~m - =  ~m -  (3) 

J2nk47trk  J2^rk  J2^k/rk 

k= 1  k= 1  k= 1 


where  rk  are  the  particle-radii.  The  volume  fractions  may  be  repre¬ 
sented  as 


M  M 

J2n'e47tre/3  Y,nVi 

1= 1  1= 1 


where  n£  is  the  number  of  /c-particles  per  unit  volume.  The  num¬ 
ber  fractions  represent  the  fraction  of  k-particles,  which  can  be 
written  as 


Sk  = 


fk/rl 


(5) 


l=\  1= 1 


The  volume  fractions  \j/k  are  defined  relative  to  the  total  solid  mate¬ 
rials,  with  the  gas-phase  porosity  0g  specified  independently.  For 
example  a  composite  electrode  that  is  (by  volume)  33%  Ni,  33%  YSZ, 
and  34%  pores  has  V'Ni  =  V^ysz  =  0.5. 

The  average  coordination  number  of  all  particles  within  an  M- 
component  packed  bed  can  be  expressed  as 


M 

z = 

k=  1 


(6) 


where  Zk  represents  the  average  contact  number  between  fc-type- 
particles  and  all  neighboring  particles  (Eq.  (1)). 

The  variable  Nk  i  (Eq.  (2))  is  particular  to  a  certain  idealized 
two-component  particle  mixture  (/<  and  l  particles)  [18].  This  ide¬ 
alization  rests  upon  two  assumptions.  First,  the  mixture  is  assumed 
to  be  dominated  by  i  particles  with  very  few  /c-particles.  In  other 
words,  ijk  (the  fraction  of  fc-particles)  is  nearly  zero  and  &  (the 
fraction  of  i  particles)  is  nearly  unity.  Second,  the  overall  average 
coordination  number  within  the  idealized  medium  is  the  same  as  it 
is  within  an  M-component  mixture.  With  these  assumptions,  Nki 
is  a  function  of  only  Z  and  the  ratio  of  the  particle-radii.  In  the  case 
that  all  particles  are  the  same  (i.e.,  k  =  £),  the  relationship  between 
Nk  i  and  Z  is 


N/a=z 


(7) 


Combining  Eqs.  (2)(3)  and  (7),  the  coordination  number  among  the 
same  kind  of  particles  Zk  k  is  a  function  of  the  volume  fraction  and 
particle-radii  as 


Zk,k 


2  fkllK 

M 


i 


(8) 


This  expression  leads  to  a  good  estimate  in  binary  particle  mixtures 
of  the  relationship  between  the  percolation  threshold  for  same- 
type  particles  (i.e.,  percolation  threshold  for  k  or  l  particles)  and 
the  particle-size  ratio  (i.e.,  rfe/r^)  [12,26]. 

In  cases  where  k  =/=  l,  the  relationship  between  the  Nk i  and  Z  is 
more  complicated  than  is  indicated  by  Eq.  (7).  Suzuki  and  Oshina 
proposed  that  Nk ^  can  be  estimated  as  [18], 

„  0.5(2  -  V3)Nkj<(rk/ri  +  1) 

ivk  i  = - T-75-.  fyj 

1  +  -  (rfc/r€(rk/r^  +  2))  7 


Fig.  1.  Relative  error  in  contact  conservation  associated  with  applying  Eq.  (9)  in  a 
multi-component  particle  mixture. 


However,  this  expression  does  not  satisfy  the  necessary  contact- 
number  conservation  principle  (i.e.,  f kZki  =  &Zgjk).  Contact- 
number  conservation  means  that  the  net  contact  between  the 
k-particles  and  ^-particles  must  be  equal  to  the  contact  between  l- 
particles  and  k-particles  [12].  The  relative  difference  between  £kZk  j 
and  ifiZg  k  when  using  Eq.  (9)  can  be  evaluated  as 

_  I Skzk,i  ~  Sizi,k\ 

min {tJkzk,l’  Sizi,k) 

This  relative  error  is  a  function  of  the  particle-radii  ratio  (r^/rfc), 
but  is  independent  of  the  average  coordination  number  Z  and  the 
number  of  components  M  in  the  mixture.  As  illustrated  in  Fig.  1, 
the  relative  error  increases  greatly  as  the  particle  sizes  vary.  Thus, 
an  alternative  approach  is  needed  to  determine  Nk 

Beginning  with  some  fundamental  principles  of  percolation  the¬ 
ory,  an  alternative  approach  is  taken  for  evaluating  Nk^.  For  a  binary 
mixture,  the  following  set  of  equations  must  hold: 


Zk  =Zk,k  +Zk,h  Zt=Zl,tJrZl,ki  (H) 

Z  =  SkZk  +  SlZl,  (12) 

SkZk,i  =  SiZl,k>  (13) 

ft  +  ft  =  l,  (14) 


Zk,k  —  Nj 


%k 


Sk rl 


m+tin 


zi,i = 


sri 


sri + sri 


Zk,l  =  SlNk,l 


(15) 


(16) 


Substituting  Eqs.  (11)  and  (13)  into  Eq.  (12),  the  average  coordina¬ 
tion  number  can  be  written  as 


Z  =  SkZk,k  +  StZt,t  +  2SkZk,l- 


(17) 


The  coordination  number  Zk  t  can  be  determined  by  substituting 
Eq.  (15)  into  Eq.  (17),  yielding 


=  0.5(1  +rl/r2l)StZ 


z  z  (, 

By  comparing  Eq.  (18)  with  Eq.  (16),  it  is  evident  that 
NW  =  0.5(1  +r2k/rj)Z. 


(18) 


(19) 


Eq.  (19)  is  equally  applicable  for  k  =  i  and  k  /  L  It  is  clear  that 
Nk  i  is  a  function  of  only  Z  and  the  radii  of  the  k-  and  ^-particles. 
Consequently,  Eq.  (19)  can  be  used  to  evaluate  the  coordination 
number  between  two  particle  types  in  a  multi-component  mixture. 

Using  Eq.  (4),  Eq.  (18)  can  be  rewritten  in  terms  of  the  particle 
volume  fractions  ij/k  as 


ZM  =  0.5(1  +r2k/rj)Z^^  .  (20) 
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Electrode  interconnect 


Dense  electrolyte 

Pel  — 

\  / 3.764  -Zel,el\2  5' 

V  2  ) 

V 

_ 2 

(i.eM  from  the  dense  electrolyte  to  the  electrode  current  collector). 
B  clusters:  The  electrolyte  particles  form  a  short  network  that  is 
connected  only  to  the  dense  electrolyte.  Such  clusters  can  carry 
ionic  current  into  the  composite-electrode  structure,  assuming 
that  there  is  sufficient  connectivity  within  the  electrode  parti¬ 
cles  to  carry  electronic  current  and  thus  facilitate  charge-transfer 
reactions. 

C  clusters:  The  electrolyte  particles  form  a  completely  isolated 
cluster.  These  clusters  are  essentially  ineffective  in  supporting 
charge-transfer  or  current  conduction. 

The  probability  that  electrolyte  particles  belong  to  an  ^-cluster 
is  defined  as  the  percolation  probability  of  the  electrolyte  particle. 
The  percolation  probability  can  be  estimated  as  [12] 


(23) 


Fig.  2.  Conceptual  illustration  of  a  composite  SOFC  electrode.  The  dark  particles 
represent  the  electrode  particles  and  the  light  particles  represent  the  electrolyte. 


The  coordination  number  for  the  /c-particles  can  be  obtained  by 
substituting  Eqs.  (8)  and  (20)  into  Eq.  (1),  yielding 


1V1 

Zk  =  =  \ 


1  +  ■ 


X>r 


(21) 


The  validity  of  the  alternative  model  proposed  above  is  supported 
by  comparison  with  an  expression  proposed  by  Bouvard  and  Lange 
[17].  Considering  a  binary  mixture  with  an  overall  coordination 
number  Z  =  6,  Eq.  (21 )  yields 


Zk  =  3  + 


3rfc 

tel  +  tej ' 


(22) 


This  result  is  the  same  as  that  reported  by  Bouvard  and  Lange  [17]. 
Further  validation,  as  discussed  below,  is  achieved  by  comparing 
with  the  computational  simulated  results  for  binary  mixtures  as 
reported  by  Suzuki  and  Oshima  [18]. 


As  can  be  seen  from  Eq.  (8),  the  percolation  probability  depends 
upon  the  particle  volume  fraction,  the  particle-radii,  and  the  aver¬ 
age  coordination  number  Z.  There  is  a  threshold  in  the  volume 
concentration  of  the  electrolyte  particles  \jrceV  called  the  percolation 
threshold,  below  which  the  particles  form  only  B  and  C  clusters.  For 
volume  fractions  above  the  threshold,  networks  of  A  clusters  are 
formed.  Above  the  percolation  threshold,  the  composite  electrode 
still  contains  a  few  B  and  C  clusters.  However,  Costamagna  shows 
that  above  the  threshold  most  of  the  particles  belong  to  A  clusters 
[12]. 

Combining  Eqs.  (8)  and  (23)  shows  that  the  percolation  thresh¬ 
olds  for  electrode  and  electrolyte  particles  may  be  represented  in 
terms  of  the  threshold  volume  fractions  (i/Ad  and  t/^)  and  the  ratio 
of  particle-radii  (rel/red)  as 


Zed.ed  —  Z 


Ze\,el  ~  Z 


V^Aed 


V^dAed  +  fe\/re\ 


f e,Ael 


V^edAed  +  V^iAel 


.  1.764 

1.764 


(24) 

(25) 


These  relationships  assist  in  identifying  parameters  that  are  impor¬ 
tant  for  predicting  the  performance  of  composite  electrodes.  For 
example,  if  the  particle-radii  are  fixed,  then  the  volume  loading 
fractions  should  be  established  such  that 


3.  Binary  mixture  properties 

Initially  consider  applying  the  percolation  theory  to  a  two- 
component  SOFC  porous  composite  electrode  that  is  composed 
of  ionic  and  electronic-conducting  particles  (labeled  “el”  for  elec¬ 
trolyte  such  as  YSZ  and  “ed”  for  electrode  such  as  Ni).1 

3.1  Percolation  probabilities 

Fig.  2  illustrates  a  two-component  randomly  packed  mixture  of 
the  spherical  particles,  where  three  different  cluster  types  are  rep¬ 
resented  for  the  electrolyte  particles.  Following  the  nomenclature 
of  Costamanga,  et  al.  [12],  such  clusters  may  be  characterized  as 

A  clusters:  The  electrolyte  particles  form  a  percolated  cluster  that 
extends  through  the  entire  thickness  of  the  composite  electrode 


1  Throughout  this  paper  the  word  “electrode”  has  two  different,  but  commonly 
used,  meanings.  On  one  hand,  electron-conducing  particles  are  called  electrode  par¬ 
ticles.  On  the  other  hand,  the  entire  composite  structure  is  called  the  electrode.  In 
context,  the  meanings  are  clear. 


V'ed  <  V'ed  <  0  —  V'eA  (26) 

By  satisfying  this  constraint  both  the  electrode  and  electrolyte  par¬ 
ticles  should  connect  (i.e.,  percolate)  through  the  entire  thickness 
of  the  composite  electrode. 

3.2.  Effective  three-phase  boundary  length 

The  contact  perimeter  between  overlapping  electrode  and  elec¬ 
trolyte  particles  can  be  evaluated  as  lk  i  =  2nrc ,  where  rc  is  the  neck 
radius  (Fig.  3).  The  neck  radius  depends  upon  the  smaller-particle 
radius  and  the  contact  angle  0  as  [12] 

rc  =  min(rfc,  r^)sin0  (27) 

When  the  intersection  between  electrode  and  electrolyte  particles 
is  exposed  to  reactive  gases  within  pores,  the  exposed  length  of 
the  intersection  is  the  three-phase-boundary  (TPB)  length.  Here  the 
entire  intersection  length  is  assumed  to  be  TPB  length. 

The  TPB  length  per  unit  volume  of  a  composite  electrode 
with  a  binary  mixture  is  assumed  to  be  the  product  of  the  con¬ 
tact  perimeter  between  electrode  and  electrolyte  particles  and  the 
number  of  contact  points  per  unit  volume,  which  is  the  product  of 
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Fig.  3.  Overlapping  k-  and  ^-particles,  forming  a  three-phase  boundary  around  the 
intersection. 


the  particle  number  per  unit  volume  and  the  coordination  number 
[12].  That  is 

^TPB  =  ^ed.el^ed^ed.el  —  ^ed.el^el^el.ed’  (28) 

where  the  number  of  k  particles  per  unit  volume  within  the  entire 
composite-electrode  structure  can  be  estimated  as  [3,12,18] 


v  W-4>g)&k 
k  4x^/3 


(29) 


The  porosity  (i.e.,  gas-phase  volume  fraction)  is  represented  as 
(fig.  Generally  speaking,  electrodes  are  designed  to  have  large  TPB 
length.  However,  TPB  length  is  only  effective  electrochemically  if 
the  participating  particles  are  percolated  through  the  thickness  of 
electrode  structure  over  which  charge-transfer  chemistry  proceeds. 
For  typical  SOFC  electrodes,  this  charge-transfer  length  is  typically 
a  few  tens  of  microns  [11  ]. 

To  illustrate  the  effects  of  relative  particle  loading,  Fig.  4  shows 
three  different  electrode-particle-loading  scenarios  as  explained 
below: 


(a)  Low  electrode-particle  loading:  Low  electrode-particle  loading 
occurs  when  0  <  t/red  <  x//ce6.  As  illustrated  in  Fig  4  a,  most  of  the 
electrolyte  particles  are  in  A-type  clusters,  with  relatively  few  in 
B-  or  C-type  clusters.  In  this  case,  the  effective  TPB  regions  (noted 
with  the  bracket  on  the  left)  are  found  only  near  the  electrode¬ 
interconnect  surface.  These  TPBs  are  potentially  electrochemically 
active.  However,  the  practical  electrochemical  activity  depends 
upon  the  ability  of  ions  to  be  transported  from  the  dense  elec¬ 
trolyte  to  the  TPBs.  This,  in  turn,  depends  upon  the  thickness 
of  the  composite  electrode  and  the  ion  conductivity  of  the  elec¬ 
trolyte  particles.  For  very  thin  composite  structures,  where  ions 
can  propagate  through  the  entire  electrode  with  relatively  low 
potential  losses,  most  of  the  TPB  lengths  in  Fig  4  a  are  electrochem¬ 
ically  active.  Such  behavior  is  observed  experimentally,  where 
the  polarization  resistance  does  not  depend  strongly  upon  the 
composite-electrode  loading  fractions  of  electrode  and  electrolyte 
particles  [12,27].  For  relatively  thick  composite  electrodes,  ions 
cannot  propagate  significant  distances  from  the  dense-electrolyte 
surface  because  of  low  ion  conductivity  and  thus  large  ohmic 
losses.  In  this  case,  the  TPB  lengths  near  the  electrode  current 
collector  are  essentially  ineffective,  rendering  the  entire  compos¬ 
ite  electrode  to  be  ineffective.  Again,  this  behavior  is  observed 
experimentally.  In  relatively  thick  composite  electrodes,  polariza¬ 
tion  resistance  increases  sharply  when  the  electrode  particles  are 
not  fully  percolated  (i.e.,  connected  networks  through  entire  struc¬ 
ture)  [12,28,29].  Generally  speaking,  the  ohmic  resistance  (ohmic 
overpotential)  is  high  in  composite  electrodes  that  use  electrolyte 
materials  with  relatively  low  ion  conductivity  or  do  not  have  fully 
percolated  electrode  particles. 


(b)  Equivalent  particle  loading  Equivalent  particle  loading  occurs 
when  <  0ed  <  (1  -  xjrceX).  As  shown  in  Fig.  4  b,  most  of  the 
electrode  and  electrolyte  particles  are  found  in  A-type  clusters 
[12].  The  TPBs  are  distributed  throughout  the  entire  composite- 
electrode  structure.  Electrochemical  charge  transfer  can  proceed 
via  ample  TPBs  near  the  electrolyte  surface.  Although  there  is  avail¬ 
able  TPB  throughout  the  electrode  structure,  the  charge-transfer 


Dense  electrolyte 
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(b)  Electrode  interconnect 


Dense  electrolyte 
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Fig.  4.  Illustration  of  composite-electrode  structures  for  three  different  loadings 
of  electrode  (dark)  and  electrolyte  (light)  particles:  (a)  low  electrode-particle  load¬ 
ing,  0  <  i/fed  <  ^ld- (b)  Equivalent  electrode-  and  electrolyte-particle  loading,  < 
\fsed  <  (1  —  V^ei)*  (c)  electrode-particle  loading,  ( 1  -  )  c  i/red  <  1. 
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region  is  usually  concentrated  near  the  dense  electrolyte  (typi¬ 
cally  tens  of  microns).  Because  of  relatively  low  ion  conductivity 
in  the  electrolyte  particles  (e.g.,  YSZ),  the  available  charge  is  trans¬ 
ferred  from  O2- within  the  YSZ  to  electrons  in  the  conduction  band 
of  the  electrode  particles  (e.g.,  Ni)  over  a  relatively  short  distance 
(i.e.,  the  electrochemically  active  zone).  Beyond  the  electrochemi- 
cally  active  zone,  the  electrolyte  particle  clustering  is  unimportant. 
However,  it  is  important  the  electrode  particles  percolate  fully  to 
carry  electric  current  up  to  the  electrode  interconnect.  The  results 
in  this  paper  are  concerned  primarily  with  composite  electrodes 
having  equivalent  particle  loading. 

(c)  High  electrode  particle  loading:  High  electrode-particle  loading 
occurs  when  (1  -  \jrceX)  <  0e d  <  1.  As  the  electrode-particle  load¬ 
ing  increases  (Fig.  4c),  most  electrode-particles  belong  to  A-type 
clusters,  with  the  electrolyte  particles  being  dominantly  in  B-  and 
C-type  clusters.  In  this  situation,  the  effective  TPBs  are  concen¬ 
trated  near  the  dense-electrolyte  surface.  This  situation  could  be 
acceptable  in  an  SOFC  anode.  The  charge  transfer  proceeds  in  the 
electrochemically  active  zone  near  the  dense  electrolyte,  and  the 
electronic  current  is  easily  conducted  upward  toward  the  anode 
interconnect.  However,  high  electrode-particle  loading  may  come 
at  the  expense  of  too  few  electrolyte  particles  and  consequently 
too  small  TPB  length.  Such  a  structure  would  cause  high  activation 
overpotentials. 

The  effective  TPB  length  per  unit  volume,  which  is  needed 
to  facilitate  electrochemical  charge  transfer,  depends  upon  the 
geometric  TPB  length  (i.e.,  independent  of  cluster  type)  and  the 
percolation  probabilities  of  the  electrode  and  electrolyte  particles. 
The  effective  TPB  length  per  unit  volume  can  be  written  as 

^TPB.eff  =  ^TPB  ^ed  ^el  •  (30) 

3.3.  Contact  area 

The  overall  contact  area  per  unit  volume  A^d  d  is  assumed  to  be 
proportional  to  the  cross-sectional  area  per  contact  between  elec¬ 
trode  and  electrolyte  particles  and  the  number  of  contact  points  per 
unit  volume.  That  is, 

^ed,el  =  ^ed,el^ed^ed,eb  (31) 

where  ak ^  =  7rr2  is  the  cross-sectional  contact  area  between  a  k- 
particle  and  a  ^-particle  (Fig.  3). 

3.4.  Effective  TPB  length  per  electrolyte  surface  area 

Electrochemical  charge  transfer  primarily  takes  place  at  three- 
phase  boundaries  that  are  distributed  within  the  electrochemically 
active  zone  near  the  dense  electrolyte.  However,  there  is  also  charge 
transfer  at  the  interface  between  the  dense-electrolyte  surface  and 
intersecting  electrode  particles  [3  ].  The  effective  TPB  length  per  unit 
electrolyte  surface  area  is  assumed  to  be  proportional  to  the  contact 
perimeter  between  an  electrode  particle  and  the  electrolyte  surface. 
This  TPB  length  may  be  expressed  as  the  product  of  the  number  of 
contact  points  per  unit  electrolyte  surface  area  and  the  percolation 
probability  of  the  electrode  particle 

^TPB.eff  =  ^ed.ele^ed^ed*  (32) 


3.5.  Electronic  and  ionic  conductivity 

It  is  generally  understood  that  there  are  important  relationships 
among  the  effective  electric  conductivity  of  the  composite  elec¬ 
trode  and  the  electrode’s  microstructure,  including  composition, 
porosity,  particle-size  ratio,  etc.  The  effective  conductivity  through 
the  electrode  can  be  separated  into  intra-particle  conductivity  and 
inter-particle  contributions  [14,15].  The  net  conductivity  is  deter¬ 
mined  as  the  result  of  effective  resistances  adding  in  series. 

3.5.1  Effective  intra-particle  conductivity 

Chen  et  al.  report  that  the  effective  intra-particle  ionic  con¬ 
ductivity  <7tra,eff  (i.e.,  conductivity  within  particles)  is  primarily 
determined  by  the  effective  relative  density  of  electrolyte  particles 
$ei  [15].  Furthermore  Jeon  et  al.  report  that  intra-particle  conduc¬ 
tivity  can  be  estimated  as  [3] 

^tra.eff 

=%  =  [(!  -faWelPeir  (34) 

°el 

where  <rda,eff  is  the  effective  intra-particle  conductivity  based  on 
the  geometric  dimensions  and  <7da,°  is  the  intrinsic  material  con¬ 
ductivity.  The  Bruggeman  factor  ju  is  used  to  include  the  effects  of 
tortuous  conduction  paths  (/z  is  typically  about  1.5).  An  expression 
analogous  to  Eq.  (34)  is  obtained  for  the  electronic  conductivity 
through  the  electrode  particles. 


3.5.2.  Effective  inter-particle  conductivity 

The  properties  at  particle  interfaces  (or  grain  boundaries)  can 
play  an  important  role  on  overall  conductivity.  Especially  at  inter¬ 
mediate  and  low  temperatures,  relatively  low  conductivity  of 
ceramic  ion  conductors  (e.g.,  YSZ)  can  be  attributed  to  low  con¬ 
ductivity  at  particle  interfaces  [15,30].  Several  factors,  including 
particle  size,  inter-particle  area,  doping  level,  and  inter-particle 
material  thickness,  can  significantly  affect  inter-particle  conduc¬ 
tivity.  According  to  Chen  et  al.  [14],  the  ion-transport  resistance  can 
be  expressed  in  terms  of  the  effective  inter-particle  conductivity 
<7dr,eff  and  the  intrinsic  conductivity  of  the  particle  interface  <7dr’°, 


nter 

Ke\ 


L 


ater,effS 

el  3 


ter,o  rter  ' 
el  Jel 


(35) 


In  this  expression  S  and  L  represent  the  geometric  cross-sectional 
area  and  thickness  of  the  composite  electrode.  Ldr  is  the  total 
thickness  of  the  particle  interfaces  along  the  current-conduction 
path.  Assuming  that  a  composite-electrode  can  be  represented 
schematically  as  in  Fig.  2,  it  is  reasonable  to  assume  that  structure 
is  composed  of  particle  layers.  In  this  case,  the  net  length  of  the 
interface  material  is  the  product  of  the  number  of  layers  and  the 
thickness  of  an  individual  particle  interface  8. 

The  inter-particle  conductivity  can  be  evaluated  as 


_ter,eff 


_ter,o 


2rdS^ 

8S 


(36) 


where  S^r  is  the  effective  contact  surface  area  between  electrolyte 
particles  per  particle  layer, 


=  aei,ei(2reiS)nV?^pel 


(37) 


In  this  expression,  “ele”  represents  the  dense  electrolyte  surface, 
and  lk  ele  =  2nrk  sin  6  is  the  contact  perimeter  between  a  k  parti¬ 
cle  and  the  electrolyte  surface.  The  number  of  k  particles  per  unit 
electrolyte  surface  area  can  be  estimated  as  [3] 


s  =  U  -<ftg  Wk 
k  27H-2/3 


(33) 


In  other  words,  Sdr  is  proportional  to  the  surface  area  per  contact 
between  two  electrolyte  particles  ael  ei,  the  number  of  the  elec¬ 
trolyte  particles  within  a  particle  layer  (2rel5)nd,  the  coordination 
number  Zelel,  and  the  percolation  probability  of  the  electrolyte  par¬ 
ticles. 

For  the  ceramic  materials,  especially  for  the  semiconducting 
materials,  the  inter-particle  conductivity  can  significantly  influence 
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the  overall  conductivity.  However,  for  the  electron-conducting  met¬ 
als  the  resistance  caused  by  the  inter-particle  is  usually  negligibly 
small. 

3.6.  Mean  hydraulic  pore  radius 


For  such  a  mixture,  the  total  volume  fractions  of  electrode  and 
electrolyte  particles  may  be  evaluated  as 

m  n 

fed  =  ^2  ^edk  ^el  =  ^2 f*'k  ■  (41 ) 

k=  1  k= 1 


The  hydraulic  radius  is  used  commonly  for  describing  flow  in 
non-circular  tubes  and  channels.  Using  a  hydraulic  radius  enables 
flow  in  non-circular  channels  to  be  analyzed  analogously  with  rel¬ 
atively  simpler  circular  channels.  The  average  pore  radius  is  also  an 
essential  parameter  in  modeling  multi-component  porous-media 
gas-transport  using  a  Dusty-Gas  model  [11,10]. 

The  hydraulic  radius  for  gas  flow  within  a  porous  electrode  can 
be  written  as 


The  relative  volume  fractions  of  particle  sizes  within  the  electrode 
and  electrolyte  materials  are  represented  as 


^eV 


f± 

fed  ’  elk  fe  1  ' 


(42) 


The  average  coordination  number  for  all  n  electrolyte  particles  is 
proposed  to  be 


2</>g 


rs  =  w  ’ 

^s.g 


(38) 


where  A^g  is  the  interfacial  area  per  unit  volume  between  the  solid 
particles  and  pores.  The  interfacial  area  can  be  estimated  as 


AlS  =  x4?r(nedred  +  nelrel )  (39) 

where  x  is  an  adjustable  factor  that  represents  the  fraction  of  solid 
surface  areas  belonging  to  the  solid-pore  binary  phase  boundary 
areas.  It  is  usually  reasonable  to  assume  that  factor  x  is  the  porosity 
0g.  Using  Eqs.  (38)  and  (39)  the  hydraulic  pore  radius  maybe  written 
as 
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fed  ^ed  +  fe\/re\ 


(40) 


n  n 

7  k=  1  1= 1 

Zel.el  =  - „ - 


k= 1 


(43) 


where  fel  is  the  number  fraction  of  e\k  particles.  The  percolation 
probability  Pel  for  electrolyte  particles  is  determined  by  the  coor¬ 
dination  numbers  among  the  same  kinds  of  particles  Zeifc?eifc  and  by 
the  mutual  coordination  numbers  Zeik?eir  U<  +  t).  Substituting  Eq. 
(43)  into  Eq.  (23)  yields  the  percolation  probability  Pe \. 

As  in  the  binary  situation,  there  is  also  a  threshold  coordina¬ 
tion  number  for  the  electrode  and  electrolyte  particles,  Zeded  and 
Zelel,  respectively.  Above  the  threshold  the  particles  are  connected 
through  the  entire  composite-electrode  thickness  [12]. 


4.  Multi-component  mixture  properties 


Zeded  =  1.764,  Ze\e\  =  1.764 


(44) 


Percolation  theory  for  multi-component  mixtures  is  more  com¬ 
plex  than  it  is  for  binary  mixtures.  This  is  especially  the  case  when 
there  is  a  distribution  of  particle-sizes.  Fig.  5  illustrates  a  composite 
electrode  with  a  single  electrode  particle  size  and  two  sizes  of  elec¬ 
trolyte  particles  (designated  eli  and  el2 ).  All  the  effective  properties 
of  the  composite  structure  depend  upon  the  particle-size  distribu¬ 
tion.  The  percolation  probability  depends  upon  the  coordination 
numbers  among  the  similar  particle  types  (i.e.,  Zell  ^ ,  Zelz  ei2 )  and 
also  the  mutual  coordination  numbers  (i.e.,  Zell  et2  and  Zel2?dl ). 

Consider  a  composite  electrode  that  includes  m  particle  sizes 
for  electronic  conductive  material  (designated  edi ,  ed2-  •  edm)  and 
n  particle  sizes  for  electrolyte  material  (designated  eli,  el2-  •  eln). 


Fig.  5.  Illustration  of  a  composite  electrode  with  two  different  electrolyte  particle 
sizes. 


Variables  contributing  to  the  percolation  threshold  are  the  radii 
ratios  and  particle  volume  fractions  i/s°dj{kem),  and 

i/red.  Assuming  that  the  radii  ratios  among  all  particle  types  and 
all  volume  fractions  corresponding  to  the  materials  (^d  anc*  V^i, ) 
are  fixed,  thresholds  in  the  total  volume  fractions  of  electrode 
and  electrolyte  particles  (^d  and  t/^)  can  be  determined.  Above 
these  thresholds  electrode  and  electrolyte  materials  are  connected 
through  the  composite-electrode  structure. 

The  effective  per  unit  volume  TPB  lengths  can  be  evaluated  as 

m  n 

ZTPB,eff  =  *edfc ,  el€  nedkZed  k ,  e\t  PedPe\ >  (45) 

k= 1  1= 1 

where  is  the  number  per  unit  volume  of  electrode  particles  in 
size  class  k.  The  contact  perimeter  between  overlapping  electrode 
(edfe)  and  electrolyte  particles  (elg)  kdk,e\£  =  27rrc  is  evaluated  as  in 
the  binary  case,  but  now  accommodating  a  range  of  particle  sizes 
(Fig.  3). 

The  effective  TPB  lengths  per  unit  electrolyte  surface  area  can 
be  estimated  as 

m 

ZTPB,eff  =  y^edk,elened,,/ech  (46) 

k= 1 


where  n is  the  number  per  unit  surface  area  of  electrode  particles 
in  size  class  k. 

The  effective  intra-particle  ion  conductivity  can  be  estimated  as 


^el 

This  expression  is  analogous  to  Eq.  (34)  for  the  binary  situation. 
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The  effective  conduction  contact  areas  between  e\k  particles  and 
neighboring  electrolyte  particles  can  be  estimated  as 

n 

SX  =  Z>Wrel  kS)n^Zelk,e,/el  (48) 

1=1 

The  effective  inter-particle  ionic  conductivity  can  be  estimated  as 


y^2relkSelk 

^.ter.eff  _  ^ter.o  k= 1 
el  -  °el  8S 


(49) 


The  mean  hydraulic  radius  in  multi-component  mixture  may  be 
estimated  as 
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(50) 


5.  Results  and  discussion 

5.1  Model  validation 

It  is  generally  difficult  to  measure  the  coordination  number  in 
either  binary  or  multi-component  mixtures.  Simulations  of  random 
arrays  of  particles  are  usually  used  to  establish  coordination  num¬ 
bers  for  particular  packing  scenarios.  Fig.  6  shows  electrode  and 
electrolyte  coordination  numbers  in  binary  mixtures  as  functions  of 
the  fractional  loading  of  electrode  particles.  In  all  cases,  Z  =  6.  Four 
different  electrolyte-electrode  particle-radii  ratios  are  considered, 
ranging  from  1.5  to  4.  Fig.  6  shows  very  good  agreement  between 
the  percolation  model  and  computations  reported  by  Suzuki  and 
Oshima  [18]. 

For  any  particular  composite-electrode  structure  Z  depends 
upon  details  of  the  fabrication  processes.  Without  further  knowl¬ 
edge  (experimental-  or  model-based)  for  a  particular  structure,  it 
is  widely  assumed  that  Z  =  6  [12,14,31,32].  As  discussed  by  Ali  et 
al.  [24],  a  model  composite  electrode  may  be  synthesized  by  first 
assuming  a  random  packing  of  spheres  followed  by  sphere  growth 
that  simulates  sintering.  The  random  packing  should  lead  to  Z  «  6, 
and  moderate  particle  growth  should  affect  the  overall  coordination 
number  only  weakly.  The  graphical  results  discussed  in  subsequent 
sections  of  the  present  paper  are  based  upon  assuming  thatZ  =  6. 

Because  the  model  is  analytic  it  can  easily  accept  different  val¬ 
ues  for  Z.  Indeed,  given  specific  knowledge  for  a  particular  class 
of  electrode  microstructures,  the  theory  can  be  applied  to  produce 
general  results.  Sensitivity  analysis  shows  that  percolation  thresh¬ 
olds  and  TPB  lengths  depend  significantly  upon  Z.  Flowever,  the 
conductivities  generally  are  weaker  functions  of  Z.  In  all  cases,  the 
qualitative  behavior  of  the  results  and  the  functional  forms  of  the 
nondimensional  properties  are  independent  of  a  particular  Z. 

5.2.  Binary  mixtures 

5.2.1.  Percolation  thresholds 

Percolation  thresholds  for  electrode  and  electrolyte  particles 
in  binary  mixtures  are  predicted  by  Eqs.  (25)  and  (24).  It  is  evi¬ 
dent  that  these  percolation  thresholds  depend  upon  the  ratio  of 
particle-radii.  Fig.  7  shows  the  relationship  between  the  perco¬ 
lation  threshold  (represented  as  electrode  volume  fractions,  ^d, 
1  —  V^ei)  an^  the  particle-radii  ratio  rel/red.  For  a  specified  rel/red, 
the  volume  fraction  of  electrode  particles  should  be  chosen  within 
the  range  {if/ced  <  <  1  -  t/^),  thus  assuring  that  both  the  elec¬ 

trode  and  electrolyte  particles  form  connected  networks  through 
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Fig.  6.  Comparison  of  coordination  numbers  predicted  by  percolation  theory  and 
direct  multi-particle  computations  reported  by  Suzuki  and  Oshima  [18]. 


the  entire  composite-electrode  thickness.  Fig.  7  shows  that  separa¬ 
tion  between  the  two  curves  is  greatest  near  a  particle-radii  ratio 
near  unity.  This  means  that  the  acceptable  range  of  electrode  vol¬ 
ume  fractions  becomes  smaller  as  the  difference  in  particle-radii 
increase.  It  is  also  evident  from  Fig.  7  that  both  curves  increase  as 


Fig.  7.  Percolation  thresholds  for  electrode  and  electrolyte  particles  as  functions  of 
the  particle-radii  ratio. 
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Fig.  8.  Nondimensional  effective  TPB  lengths  per  unit  volume  as  functions  of  elec¬ 
trode  volume  loading  yred  and  the  particle-radii  ratio  rei/red. 


Fig.  9.  Nondimensional  effective  TPB  lengths  per  unit  electrolyte  surface  area  as 
functions  of  electrode  volume  loading  \J/eA  and  the  particle-radii  ratio  rei/red. 


rei/red  decreases.  This  means  that  as  the  electrode  particles  become 
larger  relative  to  the  electrolyte  particles,  increased  electrode  vol¬ 
ume  loading  is  needed  to  assure  complete  percolation. 


5.2.2.  Effective  TPB  length 

Eq.  (30)  represents  the  TPB  length  in  a  binary  mixture.  Following 
substitutions  for  contact  perimeter  and  particle  number  loading, 
the  equation  can  be  rewritten  as 


iv 

ATPB,eff  ~~ 


3 

2 


min(red,re|) 


ed 


(1  -  0g)sin0 


V'ed-Zed.el^ed^el 


(51) 


The  effective  TPB  length  depends  uponZedel,  Ped  and  Pel,  which  are 
functions  of  and  the  particle-radii  ratio  rel/red.  The  TPB  length 
can  be  written  in  dimensionless  form  as 


jv 

ATPB,eff  ~~ 


l  V 

ATPB,eff 


(1  -<£g)sin  0/rf 


ed 


(52) 


The  nondimensional  effective  TPB  length  per  unit  volume  X^pB  ff 
is  independent  of  the  particular  particle-radii  (red  and  rel),  and  the 
contact  angle  0.  However,  it  does  depend  upon  the  particle-radii 
ratio,  the  porosity  </>g. 

Fig.  8  shows  how  X^pB  ff  depends  upon  the  electrode  volume 
loading  i/red  and  the  particle-radii  ratio  rel/red.  It  is  evident  from  the 
figure  that  when  rei/red  =  1  the  maximum  ^jPBeff  is  at  ^ed  =  0-5. 
When  rel/red  >  1,  the  electrode  volume  loading  for  the  maximum 
^TPBeff  decreases.  When  rel/red  <  1,  the  electrode  volume  load¬ 
ing  for  the  maximum  X^pB  eff  increases.  The  shapes  of  the  curves 
show  that  XypB  eff  is  only  defined  within  certain  threshold  limits 
that  depend  upon  the  particle-radii  ratio.  That  is,  there  are  frac¬ 
tional  loadings  above  and  below  the  thresholds  for  which  there  is 
no  effective  TPB  length. 

In  practical  terms,  the  actual  TPB  length  must  be  evaluated  from 
the  nondimensional  function.  In  this  case,  values  must  be  assigned 
for  the  average  contact  angle  and  the  porosity.  Typical  values  for  the 
contact  angle  are  around  0  =  15°  [3,12],  but  details  depend  upon  the 
particular  composite  structure. 

The  effective  TPB  length  AVpBeff  affects  SOFC  performance 
greatly.  The  net  rate  of  electrochemical  charge-transfer  chemistry 
(i.e.,  production  of  useful  electric  current)  depends  directly  upon 
the  active  TPB  lengths.  Thus,  increasing  effective  TPB  length  reduces 
activation  overpotential  and  improves  fuel-cell  performance. 

Although  the  functional  form  of  the  nondimensional  XTPBeff  is 
independent  of  the  particle  radius  (Eq.  (52)),  the  physical  ZTPB  eff 
actually  is  inversely  proportional  to  the  particle  sizes.  Consider  an 
example  in  which  red  =  rel  =  0.25  p,m  and  =  30%.  Using  Fig.  8,  it 
follows  that  ^jPBeff  =  6.05  x  1012  m/m3.  However,  with  red  =  rel  = 
50 nm,  the  effective  TPB  length  is  increased  greatly  to  A^pB eff  = 
1.51  x  1014m/m3.  All  other  things  being  equal,  this  reduction  in 


particle  sizes  increases  current  density  by  a  factor  of  approximately 
25.  Such  behavior  is  observed  experimentally  [5-8]. 

Although  usually  less  important  practically,  it  is  also  interesting 
to  consider  the  effective  TPB  length  per  unit  surface  area  of  the 
dense  electrolyte  at  its  interface  with  the  composite  electrode  (Eq. 
(32)).  As  with  the  TPB  length  per  unit  volume,  a  nondimensional 
formulation  is  helpful;  that  is 


X3 


is 

ATPB,eff 


TPB, eff  ' 


(1  -  0g)sin0/red  ’ 


(53) 


Fig.  9  shows  how  X3pB  eff  depends  upon  the  electrode  volume  load¬ 
ing  \[/ed  and  the  particle-radii  ratio  rel/red.  As  with  X^pBeff  there  is 
a  threshold  below  which  there  is  no  effective  TPB  length  per  unit 
surface  area.  In  fact,  comparing  Figs.  8  and  9  shows  that  the  lower 
thresholds  for  ^ed  are  the  same  in  both  cases.  However,  there  is 
no  upper  threshold  for  X3pB  eff.  For  a  given  electrode  volume  load¬ 
ing  fraction  (within  the  threshold),  X3pBeff  generally  increases  as 
rei/red  increases. 


5.2.3.  Effective  conductivity 

Fig.  10  shows  the  effective  intra-particle  electronic  and  ionic 
dimensionless  conductivities  as  functions  of  electrode  volume 
loading  t/red  and  the  particle-radii  ratio  rel/red.  As  in  the  case  of 
TPB  length,  there  is  benefit  in  representing  the  results  in  nondi- 


of 

Tb° 
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Fig.  10.  Nondimensional  effective  intra-particle  conductivity  for  binary  mixtures 
as  functions  of  electrode  volume  loading  \J/ed  and  the  particle-radii  ratio  rei/red.  The 
upper  panel  shows  the  electrode  conductivity  and  the  lower  panel  shows  electrolyte 
conductivity. 
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mensional  terms.  In  this  case,  comparison  with  Eq.  (34)  indicates 
the  dimensionless  form  as 


~tra,eff 

CTed 


_tra,eff 

^ed 


^tra.o 

^ed 


(i-0g  r 


[V^ed^edl^- 


(54) 


An  analogous  definition  is  used  for  the  ionic  conductivity  through 
the  electrolyte  particles,  a^a,eff. 

For  a  given  electrode  volume  fraction  (within  the  percolation 
threshold),  as  the  ratio  rel/red  increases,  d^a,eff  increases  (Fig.  10a). 
However,  as  rel/red  increases,  o^a,eff  decreases  (Fig.  10b).  At  a  cer¬ 
tain  electrode  volume  loading,  electrode  conductivity  increases 
as  the  electrolyte  particle  radius  relative  to  the  electrode-particle 
radius  increases  (as  long  as  the  particle  ratio  is  within  the  perco¬ 
lation  threshold).  However,  as  the  electrode  electric  conductivity 
increases  as  a  result  of  larger  electrolyte  particles,  the  ion  con¬ 
ductivity  decreases.  Because  of  the  larger  intrinsic  conductivity 
of  electrode  particles  relative  to  electrolyte  particles,  the  overall 
conductivity  of  the  composite  electrode  increases  with  the  increas¬ 
ing  electrolyte  particle-radii.  This  behavior  is  consistent  with  the 
results  reported  by  Yu  et  al.  [2]. 

Fig.  11  shows  the  dimensionless  inter-particle  electrode  and 
electrolyte  conductivities,  which  are  defined  as 


ter.eff 

oter’eif  = _ ?d _  (55) 

6d  CTedr'°red  Sin2  “  0g)/5  ’ 

where  8  is  the  thickness  of  any  inter-particle  material,  such  as  a 
grain  boundary.  For  YSZ,  grain-boundary  thicknesses  have  been 
estimated  to  be  around  <5^5  nm  [33]. 

Based  upon  earlier  literature,  Chan  et  al.  [14]  reported  specific 
intra-particle  and  inter-particle  electron  resistivities  for  LSM  at 
945  °C.  The  electron  conductivities  in  LSM  (nominally  the  electron 
conductor)  are  both  around  cr^a,°  =  a ^r,°  =  10000  S/m  (specific 
resistivity  around  0.01  £2  cm).  As  an  example  consider  a  situation 
in  which  T^ed  =  50%,  rel/red  =  1,  red  =  0.5  (Jim,  0g  =  0.26  and  6  = 
15°.  In  this  case,  based  upon  the  results  in  Figs.  10  and  11,  the 
effective  intra-particle  and  inter-particle  electrode  conductivities 
are  cr^a,eff  =  2126.5  S/m  and  cr^r,eff  =  99440  S/m.  For  this  material 
a^r,eff »  cr^a,eff,  which  means  that  the  inter-particle  conductivity 


Fig.  11.  Nondimensional  effective  inter-particle  conductivity  for  binary  mixtures 
as  functions  of  electrode  volume  loading  \J/ed  and  the  particle-radii  ratio  rei/red.  The 
upper  panel  shows  the  electrode  conductivity  and  the  lower  panel  shows  electrolyte 
conductivity. 


Fig.  12.  Nondimensional  hydraulic  pore-radii  as  functions  of  electrode  volume  load¬ 
ing  \fed  and  the  particle-radii  ratio  rei/red- 


is  sufficiently  large  that  it  is  essentially  negligible  in  offering  resis¬ 
tance  to  electron  conduction. 

Chan  et  al.  [14]  also  reported  ion  resistivities  for  YSZ  at  945  °C. 
The  specific  intra-particle  and  inter-particle  ionic  conductivities 
for  YSZ  are  (T^a,°  =  6.7  S/m  (specific  resistivity  around  15  fi-cm), 
and  or^r’°  =  0.05  S/m  (specific  resistivity  around  2000  £2cm).  Again 
consider  an  example  where  =  50%,  rel/red  =  1,  red  =  0.5  p,m, 
0g  =  0.26  and#  =  15°.  Under  these  circumstances,  <r^a,eff  =  1.4  S/m 
and  cr^a,eff  =  0.5  S/m.  In  this  case,  both  intra-  and  inter-particle  con¬ 
ductivity  contribute  significantly  to  the  net  ion  conduction. 


5.2.4.  Pore  hydraulic  radius 

The  pore  radius  can  be  an  important  parameter  in  the  control¬ 
ling  gas  transport  within  the  porous  electrode  structure.  Fig.  12  is  a 
nondimensional  representation  of  pore  radius  as  a  function  of  elec¬ 
trode  volume  fraction  and  particle-radii  ratio.  The  dimensionless 
radius  is  defined  as 


rg 

re d/U  -<t>sY 


(56) 


For  rel/red  =  1.0,  rg  is  independent  of  the  electrode  volume  frac¬ 
tion.  For  relatively  larger  electrode  particles  (i.e.,  rel/red  <  0.5, 
the  dimensionless  pore  radius  increases  as  V'ed  increases.  For 
relatively  smaller  electrode  particles  (i.e.,  rel/red  >  0.5,  the  dimen¬ 
sionless  pore  radius  decreases  as  V^ed  increases.  For  a  certain  t/'W. 
rg  increases  as  the  electrolyte-particle  radius  increases  relative  to 
the  electrode-particle  radius. 


5.3.  Three-component  mixtures 

This  section  considers  a  three-component  mixture,  comprised 
of  a  single  type  of  electron-conducting  particle  (labeled  ed)  and 
two  sizes  of  ion-conducting  particles  (labeled  eli  and  el2).  Such  an 
electrode  configuration  is  motivated  by  Itoh  et  al.,  who  found  that 
composite  Ni-YSZ  electrodes  with  both  coarse  and  fine  YSZ  particles 
prevented  the  agglomeration  of  Ni  and  improved  long-term  cell 
performance  [34]. 


5.3.1  Percolation  thresholds 

Eq.  (44)  involves  four  parameters  for  three-component  mix¬ 
tures:  red/re/l ,  rei2  /re/l ,  ^ed,  and  ).  Fig.  13  shows  the  percolation 
regions  as  functions  of  rei2/rell  and  ,  but  with  the  ratio  red/rell 
fixed  at  unity.  The  regions  between  the  upper  dashed  curve  and 
the  lower  solid  curve  (i.e.,  values  of  the  percolating  electrode¬ 
particle  volume  fraction)  are  the  regions  in  which  both  the  electrode 
and  electrolyte  particles  form  connected  networks  throughout 
the  entire  thickness  of  the  composite-electrode  structure.  High- 
performance  electrode  structures  are  usually  designed  to  be  within 
this  fully  percolated  region. 
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Fig.  13.  Percolation  thresholds  ^  and  (1  -  ^)  in  a  three-component  mixture  as 
functions  of  rei2  /re^ ,  and  with  red/rell  =  1. 

The  region  between  and  (1  -  t/^)  decreases  as  re/2  /re/l 
increases  at  fixed  volume  fraction  of  the  larger  electrolyte  particles 
.  There  is  a  relatively  weak,  but  still  significant,  influence  of  the 
volume  loading  of  the  larger  electrolyte  particles  t/^  (i.e.,  compare 
Fig.  13  a  and  13b).  Beyond  rej2/re/l  >  10  the  percolation  thresholds 
become  essentially  independent  of  the  particle-size  ratio. 
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Fig.  14.  Nondimensional  effective  TPB  lengths  per  unit  volume  in  a  three- 
component  mixture  as  functions  of  electrode  volume  fraction  \J/e d,  re\2/re^ ,  and  , 
with  red/re^  =  1. 


nondimensional  function  is  formed  as 

>  s 

rs  _  ATPB,eff 

TPB.eff-  -  0g)sin0/red 


(59) 


5.3.2.  Effective  TPB  length 

By  restricting  Eq.  (45)  to  a  three-component  mixture  (one  size 
electrode  particle  red,  and  two  sizes  of  electrolyte  particles,  rell  and 
rel2 ),  the  equation  can  be  rewritten  as 

17  3  /min(red,  re)l ) 

^TPBeff  =  -0g)fedSinflx  -  “  ^  ~Zed,eU  PedPel 

Z  V 


min(red,rei2) 


Zed.el'j  ^ecTel 


'ed 


(57) 


In  this  equation  Zk  Ped  and  Pei  are  functions  of  each  component’s 
volume  fraction  and  the  particle-radii  ratios  ty/rk.  Because  both 
sizes  of  electrolyte  particles  are  connected  to  form  the  network  (i.e., 
a  percolated  A  cluster),  it  may  be  assumed  that  Pe ^  =  Pel2.  Thus,  a 
single  probability  Pel  appears  in  Eq.  (57),  and  its  value  can  be  esti¬ 
mated  using  Eqs.  (23)  and  (43).  The  effective  TPB  per  unit  volume 
is  put  into  nondimensional  form  as 


% 


\v 

ATPB,eff 


TPB.eff  ‘ 


(1  -0g)  sin  e/rjd 


(58) 


Fig.  14  shows  the  nondimensional  effective  TPB  lengths  per 
unit  volume  as  functions  of  electrode  particle  volume  loading, 
electrolyte-particle-radii  ratio,  and  the  large-radius  electrolyte  vol¬ 
ume  loading.  In  all  cases,  the  electrode-particle  radius  is  the  same  as 
the  small  electrolyte-particle  radius  (i.e.,  red/rell  =1).  The  bound¬ 
aries  of  the  percolation  thresholds  are  evident  in  Fig.  14.  Outside 
these  thresholds,  where  the  particles  are  not  fully  percolated,  any 
TPB  length  is  ineffective.  It  is  also  evident  from  the  figure  that  as  the 
size  and  the  volume  fraction  of  larger  electrolyte  particles  increase, 
the  electrode  volume  fraction  for  the  maximum  nondimen¬ 
sional  TPB  length  decreases. 

Fig.  15  shows  nondimensional  effective  TPB  lengths  per  unit  sur¬ 
face  area  of  dense  electrolyte  in  a  three-component  mixture.  The 


The  nondimensional  surface-based  TPB  length  increases  monoton- 
ically  as  the  electrode  volume  fraction  increases.  At  relatively  low 
electrode  volume  loading,  the  surface-based  TPB  length  increases 
as  the  radius  of  the  larger  electrolyte  increases  (i.e.,  rel2/rell 
increases).  For  large  electrode  volume  loading  (i.e.,  >  0.5),  the 

surface-based  TPB  length  becomes  essentially  independent  of  the 
electrolyte  size  distribution. 


Fig.  15.  Nondimensional  effective  TPB  lengths  per  unit  surface  area  of  dense  elec¬ 
trolyte  in  a  three-component  mixture  as  functions  of  electrode  volume  fraction  if/ed, 
re\2  Ael! ,  and  ^ .  with  red/rdl  =  1. 
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Fig.  16.  Nondimensional  effective  intra-particle  electrode  conductivities  in  a  three- 
component  mixture  as  functions  of  electrode  volume  fraction  \J/ed,  rei2  /rdl ,  and  i /r°  , 
with  red/rdl  =  1. 


Fig.  18.  Nondimensional  effective  inter-particle  electrode  conductivities  in  a  three- 
component  mixture  as  functions  of  electrode  volume  fraction  \J/e d,  rei2  / rdl ,  and  , 
withred/rdl  =  1. 


5.3.3.  Effective  conductivities 

Figs.  16  and  17  show  the  nondimensional  intra-particle  elec¬ 
trode  and  electrolyte  conductivities  in  a  three-component  mixture 
as  functions  of  electrode  volume  fraction  ^ed,  rel2/rell,  and  t/^. 
In  all  cases,  the  electrode-particle  radius  is  the  same  as  the 
smaller  electrolyte-particle  radius  (i.e.,  red/rell  =1).  Qualitatively 
the  behavior  is  similar  to  situation  for  binary  mixtures.  Elec¬ 
trode  conductivity  increases  monotonically  as  the  electrode  volume 
fraction  increases.  For  low  electrode-particle  loading  the  nondi¬ 
mensional  intra-particle  electrode  conductivity  increases  with 
increasing  radius  of  the  larger  electrolyte  particles.  However, 
for  ^ed  >  0.5,  the  intra-particle  conductivity  becomes  essentially 
independent  of  the  relative  sizes  of  the  electrolyte  particles.  The 
nondimensional  electrolyte  intra-particle  conductivity  decreases 
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Fig.  17.  Nondimensional  effective  intra-particle  electrolyte  conductivities  in  a  three- 
component  mixture  as  functions  of  electrode  volume  fraction  \J/ed,  re  i2  /rdl ,  and  , 

with  red/rdl  =  1. 


monotonically  as  the  electrode  volume  fraction  increases.  For  low 
electrode  volume  fraction,  the  electrolyte  intra-particle  conductiv¬ 
ity  become  essentially  independent  of  the  electrolyte-particle  size 
ratio  (Fig.  17). 

Figs.  18  and  19  show  the  nondimensional  inter-particle  elec¬ 
trode  and  electrolyte  conductivities  in  a  three-component  mixture 
as  functions  of  electrode  volume  fraction  \j/edt  rel2/rell,  and 
In  all  cases,  the  electrode-particle  radius  is  the  same  as  the 
smaller  electrolyte-particle  radius  (i.e.,  red/rell  =1).  Qualitatively 
the  behavior  is  similar  to  the  situation  for  the  inter-particle  con¬ 
ductivity  in  binary  mixtures.  That  is,  the  inter-particle  electrode 
conductivity  increases  and  the  electrolyte  conductivity  decreases 
with  increasing  electrode  volume  loading.  Because  of  its  functional 


Fig.  19.  Nondimensional  effective  inter-particle  electrolyte  conductivities  in  a  three- 
component  mixture  as  functions  of  electrode  volume  fraction  \J/edl  rd2  / rdl ,  and  , 
with  red/rdl  =  1. 
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Fig.  20.  Nondimensional  hydraulic  pore  radius  in  a  three-component  mixture  as 
functions  of  electrode  volume  fraction  d.  te^Aet, .  and  \J/°^ ,  with  red/rei 1  =  1. 


form,  the  nondimensional  inter-particle  conductivity  crter’eff  is  inde¬ 
pendent  of  the  particle  radius  (Eq.  (55)).  However,  the  physical 
inter-particle  conductivity  <7ter,eff  (i.e.,  in  dimensional  form)  will 
increase  with  the  increasing  of  particle  sizes. 


5.3.4.  Pore  hydraulic  radius 

Fig.  20  shows  the  nondimensional  hydraulic  pore  radius  in 
a  three-component  mixture.  The  nondimensional  pore  radius  is 
defined  as 


red/U  ~<Pg )' 


(60) 


The  ratio  of  electrolyte  particle  sizes  has  a  strong  influence  upon 
the  pore  radius,  especially  at  relatively  low  electrode  volume  load¬ 
ing.  The  pore  radius  increases  significantly  as  the  larger  electrolyte 
particle  size  increases. 


6.  Summary  and  conclusions 

A  percolation  theory  model  is  developed  to  predict  effective 
properties  in  SOFC  composite  electrodes.  Such  properties,  includ¬ 
ing  ion  and  electron  conductivities,  three-phase  boundary  lengths, 
and  hydraulic  pore  radii,  are  needed  to  assist  optimizing  elec¬ 
trode  structures  and  for  inclusion  into  larger  scale  models  that 
use  effective  properties  in  continuum-based  models.  The  models 
accommodate  binary  and  multi-component  mixtures  of  ion-  and 
electron-conducting  particles.  The  effective  properties  for  an  entire 
electrode  structure  depend  upon  physical  characteristics,  includ¬ 
ing  particle  size,  packing  density,  intrinsic  material  properties, 
and  porosity.  Thus  the  models  provide  a  quantitative  connection 


between  measurable  and  controllable  characteristics  at  the  particle 
scale  with  overall  performance  at  the  composite-electrode  scale.  All 
the  results  are  presented  in  nondimensional  form,  which  increases 
generality  for  application  of  the  model. 
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